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The calculation of ideal gas thermody- 
namic properties for steam to 10,000 K 
is examined. Centrifugal effects are in- 
cluded using spectroscopic data for 
the lowest vibrational levels, with ex- 
tension to higher bending levels based 
on estimates from a bending model. 
Modifications are examined for rota- 
tional and vibrational cut-off effects. 
Uncertainties in obtaining a suitably 
regularized representation of energy 



versus bond stretching vibration in ap- 
proaching the dissociation energy re- 
gion appear relevant to the reliability 
of the extrapolation. 
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Introduction 



A calculation made in 1979 of the ideal gas ther- 
modynamic properties for water [1]' is the back- 
ground for this study on improving the 
extrapolation to higher temperatures. The proce- 
dure under consideration is to calculate the internal 
partition function and its first two temperature 
derivatives, using direct summation over vibra- 
tional levels. 

Direct summation was also used for the 1979 
work. It was used not only for vibrational levels, 
but for rotational ones as well for temperatures up 
to 230 K. In the earlier work of Friedman and Haar 
[2], sums over vibrational levels were computed ac- 
cording to a procedure similar to that of Mayer and 
Mayer [3] for diatomic molecules. Such a formula- 

About the Author: Harold W. Woolley, a physi- 
cist, retired from NBS some years ago but contin- 
ues his affiliation with the Bureau as a guest worker 
in its Thermophysics Division, Center for Chemi- 
cal Engineering. 

'Figures in brackets indicate literature references. 



tion is in principle a low temperature form, 
analogous to a power series in increasing powers of 
temperature, and as such encounters questions of 
convergence in extension to very high tempera- 
ture. 

While it is not too difficult in such an approach 
to include effects of lowest order anharmonicities 
in their first power contribution, the inclusion of 
higher power contributions of even the lowest or- 
der of anharmonicities including cross product ef- 
fects and rotational dependence is more complex. 
The extension to cover still higher order terms 
such as might satisfactorily describe the approach 
to the energy of dissociation would involve too 
many terms of cross product power type for easy 
confidence in the adequacy of their enumeration. 
There is in fact no formal end to the number of 
possible orders of the expansion, and each succeed- 
ing order would be enormously more complex than 
the one before it. 

In the presentation that follows, the problem of 
the myriads of correction terms is circumvented by 
using direct summation over the three vibrational 
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quantum numbers for this triatomic molecule. The 
problems which remain are the physical ones of 
obtaining reliable energies for the vibrational levels 
and reliable rotational parameters, including the 
way in which the rotational energy behaves in the 
extension to high rotational quantum numbers. 

1) As a first step, this last item is discussed with a 
study of an empirical bending model, including its 
implications as to bending energy. 

2) The next topic to be covered is the implied 
centrifugal thermodynamic effect based on empiri- 
cal constants fitting the spectroscopically-deter- 
mined rotational levels as represented by a Watson 
Hamiltonian [4]. This applies to eight low lying vi- 
brational levels, extending up to 5331 cm '. 

3) Following this, the problem of extending the 
representation of rotational details to high energy 
and high temperature is considered, including the 
termination of rotational levels due to lack of bond 
stability. 

4) Finally, considerations are introduced pertain- 
ing to a suitable extrapolation for energies of vibra- 
tional levels in approaching the dissociation region, 
a problem that is not regarded as solved. 

The calculation range of new estimates for ideal 
gas thermodynamic functions for water is from 200 
K to 10,000 K. It uses 



Q=2,Q t exp(-G v hc/kT), 



(1) 



with v indicating all sets of values for v,, v 2 and v 3 
giving stable vibrational states. G v is the vibrational 
energy in cm" 1 units. Q r is the rotational partition 
function for the given vibrational state. This is 
taken as 

Q,=Q,°f,\dQ, (2) 

with the usual semi-classical relation 

V °=7T 1/2 {kT/hcf n (A„BvC v r 1/2 a' 1 (3) 

where the symmetry number <r is equal to 2 for the 
H 2 molecule. A„ 2?» and C v are the principal rota- 
tional constants and dQ, is a "low temperature" 
quantum correction such as that of Stripp and 
Kirkwood [5] as used in reference [1]. The centrifu- 
gal distortion and stretching effects are here repre- 
sented by the factor f„ related to the "Wilson [6] 
centrifugal effect constant p," by / v = exp(A„), 
where h r =p v T. The treatment of centrifugal ef- 
fects for eight of the lowest vibrational levels is 
based on spectroscopic data interpreted with a 
Watson type rotational Hamiltonian. The extension 
to high bending quantum number is made on the 
basis of model estimates. 



An examination is also made as to plausible mag- 
nitudes for effects of rotational cut-off in the disso- 
ciation region. Ad hoc adjustments in the approach 
to the dissociation energy region have been made 
to preserve approximate symmetry versus v, and v 3 
quantum numbers. 

For temperatures upward of 200 K, the data and 
empirical representations of Camy-Peyret et al. [7] 
could be used for seven vibrational levels above 
the ground state for which detailed parameters for 
rotational Hamiltonians are available. These in- 
clude values for the principal rotational constants 
A r , B, and C», and for G,, the energy for the vibra- 
tional level at zero rotation. Their reported values 
for G, were used also for four other vibrational 
levels. Their G v data for the (1,1,1) state were not 
used, as its partner in resonance, (0,3,1), had not 
apparently been similarly covered. A placement es- 
timate for a (0,4,0) level based on a resonance shift 
from Benedict [8] was used to complete G, values 
for the resonating triad including the (2,0,0) and 
(1,2,0) levels. For the ground vibrational state, 
however, the slightly differing results of the more 
recent analysis by Kyro [9] were accepted in the 
later calculations. The overall course of vibrational 
energies versus vibrational quantum number was 
taken to follow an empirical data fit by Benedict 
[8], but with some adjustment in the higher energy 
regions to be consistent with other data, such as 
dissociation and heats of reaction. 



The Bending Model 

On the basis of spectroscopic data it can be in- 
ferred that vibrational bending produces large ef- 
fects on /„ and special effects on the total energy. 
The lack of extensive data and the impossibility of 
making reliable long range extrapolations of di- 
rectly fitted polynomial representations of data 
have led to the present numerical exploration based 
on a simple bending model. For this, somewhat 
crude evaluations of WKBJ integrals have been 
used, based on an approximate bending potential. 

The potential U is taken as the product of an 
empirical basic function U and an empirical cor- 
rection function U a as 



U=UAL 



(4) 



The basic function U here involves two Lorentz 
type terms, 

U ^k 1 /(q i +x) + k 2 /(q 2 +x)~k 1 /q 1 -k 1 /q 2 (5) 

with the correction function U c taken as 
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U c = (\+B,x +B 2 x 2 +B 3 x" -B 4 gx") 



(6) 



where *=((|> 2 -4> < , 2 )/<f> c 2 with 4> as the angle of 
bending of H-O-H out of a straight line as shown in 
figure 1. The even power of d> in the definition of* 
gives symmetry about the angle for full barrier 
height, at <f>=0. The subscript e refers to the equi- 
librium configuration. 




Figure 1-A schematic model of the water molecule. 



The calculations were made with bonds of fixed 
length "d," taking the atoms as point masses. The 
quantum phase integral 



v 2 +\/2=m H U2 d h 



[l-(l+cos<f>) 



Mmo/m H + 2)] m lG^U(i)] W7 d<t> (7) 
may be put in the form 

v 2 + 1/2=(2 1/2 /4it)(1 +cos 4> e y m f "" 

X[l -(1+cos 4>)/(m //n H + 2)] l/2 

X{[G-£/(/)]/£j l/2 d<f> (8) 

for which conventional wave number energy units 
are convenient. The quantities m and m H repre- 
sent masses of the respective atoms. 

In the absence of rotation, U(i) is identical with 
U of eq (4). For rotation purely about the principal 
axes there are three cases: 

(i=a) : U(a)^U+J d 2 A e (l -cos <f> e )/(l-cos <f>) 



(i = b) : U(b)=U+J b 2 B e (l+cos <f> e )/(l+cos <f>) 

(i-c):U(c)=U+J c 2 C e [l+m H (m + m»r l 

cos <f> e ]/[l +m H (m +m H )-' cos <f>]. (9) 

Evaluations were made using cm" 1 energy units, 
with 

U =40266.9[l/(1 -0.21507*)- 1] 

+ 13464.5[1/(1 +0.6432*)- 1] (10) 

and 

U c = (1-0.2467* - 0. 1 526* 2 + 0. 5678* 3 

+ 1.3708* 4 )/(l + l. 7442* 4 ). (11) 

A value of 75.73 degrees was used for 4> e in adjust- 
ing between A, and B e indications of a preliminary 
data fit of Benedict [8]. 

Energy values on effective potential curves in- 
cluding rotation were approximated by 



U(a)= £/ + 20.24/ a 7(1.0001 -cos <f>) 
U(b) = U + 17.97/ t 7(1.0001 +cos <f>) 
and 

t/(c)=/7 + 9.5/,7(l+0.0593cos4>) 



(12) 

(13) 

(14) 



where 0.0001 has been added in two of the denom- 
inators to avoid accidental division by zero. 

Values of v 2 + 1/2 were computed at 17 to 40 
suitably spaced values for energy "E" or "G" up 
into the 100,000 cm" 1 region for these three cases 
of rotation about the principal axes, for several val- 
ues of each J 2 ranging from zero to 400 in the "a " 
and "c" cases and somewhat further in the "b" in- 
stance. Four point Lagrangian interpolation was 
then used to obtain energies corresponding to in- 
teger Vi values at each of the chosen J, 2 values. 
Effective values for the principal rotational con- 
stants were estimated according to 



A v =(E Ja -E )/J a 2 
B v —(Ej b —E )/J b 
C v =(£ -£ )// r 2 . 



(15) 
(16) 

(17) 



Figures 2, 3, and 4 show results from these calcula- 
tions. Extrapolation to zero rotation appears reli- 
able for B v and C v and for A, for v 2 small. 



37 



Journal of Research of the National Bureau of Standards 



700 



600 



300 



ft 400 

A K' 



CM 



-1 



300 



200 



100 



Figure 2-Dependence of the effec- 
tive rotational constant Ak or 
Au on J J or K l . 
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Figure 3-Dependence of the effective rotational constant S Jt on 



Estimates for the centrifugal effect contribution 
associated with rotation about separate principal 
axes were obtained in the following way. In the 
case of the "B" rotation, for example, and empiri- 
cal representation for B versus J b 2 was used for 
each of various v 2 values. Numerical quadratures 
were performed for a partition function contribu- 
tion as 



(Q) s =rexp(-PBJ b *)dJ t 

J 



(18) 



for various temperatures where fi=hc/kT. This 
determined an effective B according to 



5 efr =77/(4^(e)/). 



(19) 



The corresponding contribution to the Wilson cen- 
trifugal effect constant then followed as 



ip^ = {2T)^\n{B /B c „). 



(20) 



A similar procedure was used for the "A " rotation. 
No appreciable contribution came from the "C" 
case. 

The combined contributions to the Wilson cen- 
trifugal effect constant based on the rigid bender 
model are shown in figure 5 for values of v 2 ranging 
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Figure 4-Effective values for the 
rotational constant Cx- 



by unit steps from 4 to 9 and also for 10, 15, 20, 25, 
and 35. (The sizes of symbols in the figure are not 
intended to indicate relative importance of the 
plotted values.) Corresponding curves are shown 
for the empirieal representation 



\0 i p=(b +b 1 t +b 2 t 2 )/(l +d\t +d 2 t 2 ), 



(21) 



where ? = 771000. For these, the numerical 
parameters that follow are based on combined esti- 
mates as from eq (21), fitted approximately by in- 
spection and graphical processes. With v 
representing v 2 , the constants b t and d, were taken 
as 

6 =2.2616(1 +0.242V +0.037v 2 + 0.00083 v 3 ) 

-Kl-O.O175v+aO033v 2 )+8/[l+O.8(v-7.3) 2 ] 

*,= l/{.122 + [-1.6+l 

-K0.31 +0.02v +0.0019v 2 - O.0OO028v 3 )] 2 } 

6 2 =1/[1.3 + 0.106(12 — v) 2 + (l£: — 06)(12— v) b 
+(2£:~08)(12-v) 8 ]-1/t5+14.3(8.3-v) 2 ] 

rf,= 12.6/[35. + (v-5) 2 ] 

rf 2 =:3./[25.+(v-ll) 2 ] 

as an approximate representation of the calculated 
values. 

Values for A„ B„ and C v as interpreted via data 
of figures 2, 3, and 4 are shown in figures 6, 7, and 
8 by solid circles. The solid curves are from a fit of 
spectroscopic data by Benedict [8] running up to 
v 2 =4, with v t and v 3 also extending up to 4. 



The x 's in figure 6 show individual estimates at 
the quantum number K — \, indicating that for A, 
extrapolation above v 2 = 7 encounters some imper- 
fection in the traditional representation. 

The dashed curves in figures 6, 7, and 8 are given 
respectively by 

A =27.8847^1 + {0.0895v +0.0228v 2 +(0.0022v 2 
+ 0.00012v'-0.000185v 4 +0.35£:-04v 5 
-0.1E -06v 6 )/[(l -0.17V +0.008v 2 ) 

X(l-0.15v+O.O07v 2 ) 2 +Lri-O6v 5 ]}) (22) 

B = 14.51 1 8{(1.0012771 +0.05266v +0.00779v 2 ) 

-Kl+0.04v+0.008v 2 ) + (2.51£:-08v 3 + 6.4£: 
-06v 4 +6.3£:-06v 5 )/(l +0-04v +0.008v 2 ) 

+0.009(c -7)/[(0.00032) 2 + (v -7) 2 ]} (23) 

C = [9.2806+0.0073v(v-l)]/(l+#+# 2 +# 3 ) 

(24) 

with i/=0.1473v/9.2806, where v represents v 2 . 

The expressions for A, and B, are roughly repre- 
sentative of the directly indicated model results as 
shown by the solid circles. For C„ the equation 
here is basically a rearrangement of Benedict's 
equation, although the model results suggest that a 
different eurve might be better. 

The open circles show values that have been ad- 
justed from the solid circles according to an al- 
lowance for bond stretching based on OH bond 
data. The effects may be summarized in part as due 
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T/110G 



Figure 5-Wilson stretching con- 
stant for H 2 versus tempera- 
ture and bending vibrational 
quantum number. Curves are by 
eq (21), points by quadratures. 
Values of v 2 for points are 



4, 

O 

5, 

6, 

o 

7, 
8, 



o 

9, 

10, 

o 

15, 



20, 



25, 
35, 



to the "effective" v 2 value, v 2 (eff), being less than 
the true v 2 in accord with 

v 2 (efO = v 2 /[l+(1.60£-03)v 2 + (6.12£-05)v 2 2 

-(5.815 -07)v 2 3 + (9.7£ -09)v 2 3 ]. (25) 

This causes the adjusted value for the "true" v cor- 
responding to each model estimate to occur for a 
larger value than the v (eff). The change in bond 
length also reduces the magnitudes of principal ro- 



tational constants from the model according to the 
square of the same ratio. 

An indication via the bending model as to the 
dependence of vibrational energy on the v 2 quan- 
tum number is also given directly by the v 2 +l/2 
values versus E values from the WKBJ integrals 
with /, 2 = 0. The results are shown graphically in 
figure 9 by the large open circles, obtained directly 
using the bending model without any bond stretch- 
ing allowance. Ad hoc adjustments for bond 



40 



Journal of Research of the National Bureau of Standards 



nn 



810 



600 



CM" 1 



too ■ 



100 ■ 




Figure 6-Estimates for the rota- 
tional constant A, 2 versus v 2 — 
solid circles via rigid bonds, 
open circles with bond stretch. 
The dashed curve is by eq (22), 
the solid curve via Benedict. 



20 



stretching were then taken from v 2 +\/2 differ- 
ences between the small open circles for free rota- 
tion without bond stretch and the small solid 
circles for free rotation with bond stretch obtained 
with bond data transferred from the OH bond of 
the OH diatomic molecule. 

As to the curves of figure 9, the one labelled 
"a" is for three terms in v 2 when v, = v 3 = 0, from a 
fit by Benedict with 2v,'s up to 4. Curve "c," as 

G = 1608.034v 2 - 1 1.748v, 2 - 1.643v 2 3 + 0.0937v 2 4 , 

(26) 



is obtained versus v 2 alone from basically the same 
data with a fourth term included in the fitting. 
Curve "b," used in the ideal gas calculations of 
1979 [1], was obtained from "a" by adding the two 
terms 

0.05v 2 (vj-l)(v 2 -2)(v 2 -3) 

- 0.0005 1 vj(v 2 - 1)0, - 2)(v, - 3)(v 2 - 4). 

Curve "d" is represented by a rational function 
with coefficients chosen to fit the large open circles 
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Figure 7-Rotational constant B n 
for rotation about intermediate 
axis - solid circles via rigid 
bonds, open circles with bond 
stretch. The dashed curve is by 
eq (23), the solid curve via 
Benedict. 




Figure 8-Estimates for 
the rotational constant 
C V2 versus Vj solid cir- 
cles via rigid bonds, 
open circles with bond 
stretch. The dashed 
curve is by eq (24), the 
solid curve via Bene- 
dict. 
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100,000 




10,000 



Figure 9- Vibrational bending en- 
ergy, Small circles: free rota- 
tion. Large circles: with 
bending potential. Open circles: 
rigid bonds. Solid circles: with 
bond stretch. Curves: (a) 3 term; 
(b) with 2 terms added to (a); (c) 
4 terra fit, eq (26); (d) rational 
function; (e) for (d) with stretch; 
(f) eombined locus, eq (27). 



v 2 



from the potential model without bond stretch. 
Curve "e" is also by a rational function, but fitted 
only to the large solid circles at low quantum num- 
ber. The last curve, "f," involves a combined locus 
asymptotic to a rational function curve at low v 2 
and to a straight line at large v 2 . It is given by 



G=(a+b)/2 + [(a-by+4c 2 ] in /2-d 



(27) 



where a=2300(v -7), c = 1600, d = 157.4661 132, 
and b = 1601.337v(l - 0.4105368v + 0.0706926v 2 



-0.006001528v 3 + 0.000229554oVH(l -0.4028641 
v+O.06850O98v 2 -0.0O5819278v 3 +0.OOO2301345v 4 ) 
where v represents v 2 . It is this curve that repre- 
sents values used in the present calculation of ther- 
modynamic functions. 

In comment on the many digits used for these 
constants, this is to be attributed not to any extreme 
accuracy, obviously, but at least in part to a regard 
for correlation between coefficients and to a desire 
to retain significance in difference type effects. 
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It is perhaps well to admit at this point that great 
accuracy is not claimed for the bending potential 
used. The application of the bending model is seen 
as quite successful, however, in providing a clear 
indication of the rather moderate magnitude of 
change in distant extrapolation as compared with 
the results for free rotation. Still, the uncertainty in 
these extrapolations must be very appreciable. 



Centrifugal Data 

As shown in an earlier publication [1], the effects 
on the rotational partition function for a given vi- 
brational state due to centrifugal distortion and 
bond stretching may be obtained in semi-classical 
approximation from the integral 

Q = 8?r 2 h ~ 3 J | Lxp(-H/kT)dP x dP y dP 2 (28) 

with H=hc(W +H l ), where ff =2B,/ > , 2 and 
where H, represents the remaining part of the Wat- 
son-type rotational Hamiltonian. This uses B 1 = C r , 
B 2 =B„ and B)~A„ the principal rotational con- 
stants, with Pi=P x , Pi=P y and P } =P Z . 

The present application of the method has been 
carried to the evaluation of five coefficients versus 
temperature. In the 1979 application of the method, 
the evaluation was made to three coefficients, only. 
The hope was that the added detail would provide 
a better overall representation of thermal effects 
implied by the spectroscopic data. 

There may typically be about 20 to 30 terms in 
current realizations of H u involving coefficients of 
various powers or products of powers of 
P* = P x i + p y i + p z \ "P xy i "=p x 1 -p y i , and P z 2 . The 
factor exp(—hcHi/kT) may be expanded as a Tay- 
lor series in powers of H v With 20 terms in H\, 
running generally up to the 10th power in jP,'s, but 
with a 12th power in P z , there are found to be 98 
terms in H, 2 and 35 in H? in the range through the 
12th power in P,. When the terms mp 2 and P xy 2 are 
expanded in powers of P x 2 , P/ and P 2 , the number 
of separate terms to use in evaluating Gaussian in- 
tegrals becomes quite large. For each, the inte- 
grand is a product over ; = 1, 2, and 3 of P 2 " 1 
exp(—hcB,Pi 2 /k7*). As a result, each separate term 
is of the form 



fi, = fi,° n F{n,)(kT/hcBtf 



(29) 



F(0)=1, F(l)=l/2, F(2)=3/4, F(3)=15/8, F(4) 
= 105/16, F(5) = 945/32, F(6) = 10395/64, etc. 

A computer program has been arranged for car- 
rying out the preparation of the correction factor 
as a series in powers of temperature using the em- 
pirical constants of the Watson-type Hamiltonian. 
With terms up to the 5th power of H x covered, the 
corresponding coefficients are computed to give 
the logarithm of the correction factor as a series 



h r =p T(\ +a, T + a 2 T^+a, T 3 +a A T% 



(30) 



where F{n ( ) = 2 ln < (2n,)\/n,\. Q° is for a classical 
rigid rotator for the level (v b v 2 , v 3 ), symbolized by 
"v" as indicated earlier. The F(«,) constants are 
simple fractions that are functions of «,■ such as 



In practical application, this has not appeared to 
provide a wcll-bchavcd form when used for mod- 
erately high temperatures. Accordingly, the pro- 
gram next computes the corresponding coefficients 
in a Pade approximant or rational function form, 

h,=p T(\ +C, T+c 2 T-)/(\ +d t T+d 2 T 1 ), 



h r =(b T + bt T 2 + b 2 r J )/(l+rf, T+d 2 r 2 ). (31) 

The coefficients in the Pade form follow from 
those preceding according to the relations 

d l =(a 1 a A -a 2 a i )/(a 1 1 -a l a } ) 

d2 = W-ai ai)/(ai 2 -a^ a } ) 

Ci—ai+di 

c, = a 2 + a l d l +di.. 

The Pade form appears to be much better 
adapted to computation in ordinary circumstances. 
For some higher vibrational states, however, there 
can still be complications such as the occurrence of 
negative values for the coefficients c 2 and d 2 which 
are for the highest powers of T. Happenings of this 
type appear to be somewhat dependent on the 
source of the empirical Hamiltonian constants 
used. 

The program is interactive in asking for values 
for the principal rotational constants, A , B, and C, 
and then for the highest power of temperature to 
be covered (up to 6 as set up). It then asks progres- 
sively, in a selected order for values for 30 coeffi- 
cients in the Watson Hamiltonian, identifying each 
by a coefficient name in "string-variable" form, 
e. g., DELJ, etc. These happen to be in the order 
shown schematically by: 

H= -DELJ*J**4-DELJK*JZ**2*J**2 

-DELK*JZ**4-2*SDLJ*JXY**2 
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-SDLK*(JZ**2*JXY**2 
+JXY**2*JZ**2) + HJ*J**6 
+ HJK*JZ**2**J**4 + HKJ*JZ**4*J**2 
+ HK*JZ**6 + 2*SHJ*JXY**2*J**4 
+ SHJK*(JZ**2*JXY**2 
+JXY**2*JZ**2)*J**2 
+ SHK*(JZ**4*JXY**2 + JXY**2*JZ**4) 
+CLJ*J**8 + CLK*JZ**8 
+ CLKKJ*JZ**6*J**2 
+ CLJK*JZ**4*J**4 + CLJJK*JZ**2*J**6 
+ 2*SLJ*JXY**2*J**6 
+ SLK»(JZ**6*JXY**2 + JXY**2*JZ**6) 
+ SLKJ*(JZ**4*JXY**2 
+JXY**2*JZ**4)*J**2 
+ SLJK*(JZ**2*JXY**2 
+JXY**2*JZ**2)*J**4 + CPK*JZ**10 
+ CPKKJ*JZ**8*J**2 
+ CPKJ*JZ**6*J**4 

+ SPK*(JZ**8*JXY**2+JXY**2*JZ**8) 
+ Z12*JZ**12 + Z10P2*JZ**10*J**2 
+ Z14*JZ**14 + Z16*JZ**16 
+ Z18*JZ**18, (32) 

Table 1. Constants for low-lying vibrational states from Watson-type Hamiltonian data. 



where J**2 is JX**2+JY**2 + JZ**2 and JXY**2 
is JX**2-JY**2. 

Values for the Pade constants as obtained from 
the available Watson Hamiltonian constants for the 
eight observed vibrational levels, based largely on 
the work of Camy-Peyret and Flaud [7] are given 
in table 1. 

A listing of the program is included in the ap- 
pendix. Further discussion of the results will be re- 
served for a later section dealing with table 
comparisons. 



Rotation at High Temperature 

While the Pade form for centrifugal effect seems 
better adapted for calculation than the simple 
power series form which encounters convergence 
problems of an erratically varying sign type, there 
are other consideration if extrapolation to very 
high temperature is required. 

The rotational quantum numbers can increase up 
to some limiting large values as "/„," "/*," or "J c ," 
beyond which centrifugal force would cause the 
molecule to break apart by bond rupture. The lim- 
iting rotational energy for rotation about any prin- 
cipal axis would be of a general magnitude 
indicated by 



D R =BJ t 2 



(33) 



where B t refers to A„ B y or C v , according to the 
axis involved. However, the affected moment of 
inertia at bond rupture would be appreciably in- 
creased by bond stretching over its ordinary value. 
A semi-classical OH bond model study on the prin- 
cipal rotational constants at zero vibration suggests 
that the ratio "r" between principal constants at 
maximum versus at low rotational quantum num- 
ber should be about 0. 195 as v A for A , 0.260 as r B for 
B, and 0.310 as r c for C. 



(v) 


Gv 


Av 


Br 


CV 


A0 


Al 


*2 


</l 


dl 


000 


0000.0000 


27.8806 


14.5216 


9.2777 


2.4518 


3.5538 


1.0949 


1.5266 


0.5377 


010 


1594.7450 


31.1284 


14.6875 


9.1291 


3.2111 


2.1258 


1.0101 


0.7945 


0.3886 


020 


3151.6301 


35.5867 


14.8415 


8.9745 


4.3691 


-.7474 


3.1192 


.01296 


0.7349 


100 


3657.0532 


27.1222 


14.3048 


9.1046 


2.4668 


-1.0581 


0.6651 


-.3572 


0.2392 


001 


3755.9296 


26.6480 


14.4313 


9.1382 


2.4097 


-.8892 


0.5480 


-.3130 


0.2144 


030 


4675.1750 


42.1323 


14.9714 


8.8350 


6.8534 


-4.9258 


10.668 


-.3500 


1.5033 


no 


5226.5870 


30.1712 


14.4139 


8.9520 


3.1892 


-1.0880 


1.8797 


-.2125 


0.5829 


011 


5331.2798 


29.5226 


14.6136 


8.9931 


3.0160 


-.6923 


1.1666 


-.1653 


0.4035 



This includes Pade type Wilson centrifugal effect parameters B(I) and D{I) for exp(p7"), where 

p={B(0)+T*[B(l)+T*B(2)]}/{l + T*[D(l) + T*D(2)]}, with B(0)=b0/l.0E + 5, i?(l) = M/1.0£ + 8, «(2)=A2/1.0£ + 11, D(l) 
= rfl/1.0£ + 3, Z>(2)=<*2/1.0£+6 
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This leads to the interesting inference that if a 
limiting partition function would be equal to the 
"volume" of an ellipsoid with semi-axes /„, J b and 
J c , one may estimate the volume as 

& w = (4/3)77-or ^(D R /A r ) [/2 (D R /B v y /2 

X(.D R /C f ) ]/1 R (34) 

where R , the factor of centrifugal increase, is esti- 
mated as 



R=(.TaT b Tc)- w ^.O 



(35) 



for low vibrational states. For high vibrational 
states, where little additional rotational energy is 
needed to bring about bond breaking, the ratio 
needed may be much nearer to unity. A form 



K=[1 + (K -1)(ZVA) S ] 



(36) 



with i? =8.0 and possibly 5 = 1 may be a useful 
speculation as to plausible behavior. Here D, is the 
additional energy D —G,, to reach dissociation for 
the vibrational state without rotation. 

As to acceptable values for D R for higher vibra- 
tional states, it appears useful, with D M as dissocia- 
tion energy including rotation, to note that by a 
logarithmic plot of y=(D M —D)/D versus 
x=(D—G v )/D for the OH diatomic potential, an 
approximate representation is y={2/9)x y> . Esti- 
mates of a similar magnitude can also be found 
from the expression y = (\/A)x , which is a form eas- 
ier to use. The latter choice provides an approxi- 
mate relation D R ^D(x+y), or D=D~~G, 
+ {D -G,)/4=\.25(D -G,). 

As in the discussion leading to eq (34), an esti- 
mate for the rotational partition function may be 
based on an integral using an ellipsoidal shell with 
semi-axes fiE/Af 1 , ft{E/B) m and UE/Cf 1 . The 
"volume" of the shell between energies hcE and 
hc(E + dE) is 2tt fi l {ABC) m E m dE. In the 

evaluation of density of states as measured by II 
(A -1 dpidqi), there is a factor 4n for orientation of 
the total momentum vector and 2ir for position of 
the rotator in making one revolution. Thus the 
number of states available within the energy shell 
needs the factor 87r 2 /!" 3 to be included, giving 



dN = 2(ABC) m E m dE. 



(37) 



An additional factor f=(}+1r i E t ) may represent 
the increase due to centrifugal effects. Integration 
to infinite energy gives 

QS =^dQ R =2(ABC) w j exp(-hcE/kT) 



X(l+2r,E')£ 1/2 rf£ 



(38) 



or 



Qi =2(kT/hcY / \ABCy 



■/: 



exp(— x)x 



(39) 



Xll + ZnikT/hcyx'jdx 

The result is 

Q S = ff w {kT/hcf n {ABC)- m 

X[l + (3/2)r,(A;r// I c)+(3/2)(5/2)r 2 (fcr/fe) 2 

+ (3/2)(5/2)a/2)r 3 (kT/hcy +...] (40) 

Using r, ={n[2/(2j+l)]}e,, this may be identi- 
fied with eq (41) 

Qj? =ir l/2 (kT/hcy n (ABCy 1/2 [l+le>(kT/hcy]. 

(41) 

The last factor is /» or exp(pT) for centrifugal 
effects according to eq (2). 

If the integration is extended only to a rotational 
energy E=D, or for x to x=hcD/kT, the result 
for each term involves an incomplete gamma func- 
tion. 

Qg = 2(kT/hc) in {ABC) " 1/2 [y(3/2,x,) 

+ 2r,y(i + 3/2,xO(kT/hcy]. (42) 

The recurrence relation y(a + l,x)=ay(ajc) 
— x a exp(— x) is used to relate all later terms to the 
first one. Values can be found for -y(3/2,jc) 
= (n l/2 /2) H(z)— z exp(— z 2 ), where x=z 2 , using 
H(z)=l-77--' /2 exp(-z 2 )(z+0.2)zl 4 [l-zl 4 (2.5803 
-2.8136 zl+4.0745 zl 2 — 1.2142 zl 3 +1.1657 zl 4 
-0.0091 z\% where zl = l/(z+0.1) ,/2 . 
The result for a given rotational state may be writ- 
ten as 



Qr d =Q\QJ,-Qs) 



(43) 



with Q°=n m (kT/hcy /2 {A v B,C r y m cr" 1 , as for a 
rigid rotator, or being the symmetry number, with 

e, = 27r- 1/2 y(3/2,x,) 

and with/, representing exp(p7") as in eq (2) or as 
used for Qg. Thus/ might be used in any form that 
would appear suitable, such as with a Pade approx- 
imant, if acceptable. Q s is a residual quantity 
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Q s = 2ir- V1 x x m exp(-jc,) le.ikT/hc)' 

x(i{n[2/(2/ + i)]}*A 

a form capable of further examination. 

As a variant study based on the derivation lead- 
ing to eq (38), one may remove the Boltzmann fac- 
tor exp(—hcE/kT) and consider the integration 
up to an energy E. 



Q E =2{A v B,C,y 



! (l+le^E 1 
Jo 



'dE. (44) 



This was tried on the ground state and some others. 
Conversion of the polynomial to a Pade- Wilson ex- 
ponential form 



f=exp[(bOE + bl E 2 + b2E } ) 
-r{\+dlE+d2E 2 )} 



(45) 



quieted a term-wise sign fluctuation effect. How- 
ever, in extension to very large E, a condition of 
excessively large computed Q E was encountered. 
This was due to the exponential factor becoming 
grossly over-sized. For the ground state, as E rises 
from l.E + 5 to l.E + 6 cm" 1 , the computed cen- 
trifugal factor rises from 6.5 to over l.E + 1. This 
result is contrary to the previous estimate of a limit 
for the centrifugal factor of the order of R = 8 or 
less, as in eq (36). The catastrophe can obviously be 
avoided by enlarging the denominator by including 
a term c/3 E i , with the ratio bl/dl near to 2 or to 
In/?. 

Logically, the parameters should be chosen 
again in such a way that the expansion into a power 
series would remain unchanged through the first 
five coefficients. If s represents b2/d3, and with 
B 1, B2, D 1 and D2 to represent original values of 
bl, b2,dl and d2, respectively, the revised coeffi- 
cients can be obtained from 

5:1=61**3 2*b0*bl*b2\b0**2*b2*dl 
-b0*bl**2*d 1 +b0**2*b l*d2 

s2=b0*b\*d\*d2-bl**2*b2-b0*b2*dl**2 
-b\*b2*d\+2*b0*h2*d2-b2**2 
-b0**2*d2**2 

b2=B2/[l+(sl/s2)/s] 

dl^bl/s 



d2 = b2*{D2-b\/$)/B2 

d\ = [bl*{D\-bQ/s)-dl*B\ 

+D2*(Bl-bO*D\)]/(B2-bO*D2) 
b\=Bl-bO*Dl+bO*dl 

With parameters so modified, the computed cen- 
trifugal factor for the ground state at l.E + 5 and 
l.£ + 6 reciprocal centimeters showed reductions 
to 2.73 and 6.54, respectively. 

The same type of control adjustment should ap- 
parently be applicable to the Pade form in terms of 
temperature in a normal computation. However, in 
actual application to a multitude of levels, there 
could seem to be a possibility that the final Pade 
constants might not always be positive, due to nu- 
merical accident. A requirement that </3 and 62 be 
positive can be met by using the absolute value of 
B2 for b2, with dl> — b2/s. The other parameters 
follow from 

bO^BO 

bl=[(Bl**2-~B0*B2)*(B\-B0*Dl) 
-B0*Bl*(B2-b2-B0*D2) 
+B0**2*(B0*d3-Dl*b2)]/DEN 

dl = [(Bl*Dl-B0*D2)*(Bl-B0*Dl) 

-B l*(B2-b2-B0*D2)+B0*(B0*d3 
-Dl*b2)]/DEN 

d2 = [(Bl*D2-B2*Dl)*(Bl-B0*Dl) 

+(B2-B0*D2)*(B2-b2-B0*D2) 

+ (B0*D\-Bl)*(B0*d3-D\*b2)]/DEN 

with DEN=(Bl**2 \ B0**2*D2-B0*B2-B0* 
Bl*Dl). 

This might preserve only four instead of five co- 
efficients of the series leading to a Pade develop- 
ment and no absolute guarantee is known to exist 
against occurrence of a zero denominator. 

Another simple scheme for keeping the rota- 
tional Q below the limiting Q M value for a given 
vibrational level has been patterned after the famil- 
iar relation of a hyperbola to its asymptotes repre- 
sented as a combination of loci. With Q u as an 
excessive rotational partition function without cut- 
off, an estimate with cut-off included might be 
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G =0/2X0* + e,)-<i/2)[4*&G A 



+ (Q M -QrfV 



(46) 



The quantity q is to be taken in a convenient form 
showing an acceptable temperature dependence. 
Results of a graphical study for rotations about 
principal axes, using bond stretching of the OH 
molecule, lead to a provisional suggestion that a 
usable form might be 



q=q x ztvp(-z) 



(47) 



where z—q 2 hcD/kT, with #,=0.4 and qi=\.%. 
Other representation schemes may reasonably be 
more suitable, however. 

It is conceded that direct rotational cut-off ef- 
fects are fairly small even for temperatures at the 
top of the range of the present tabulation. How- 
ever, an indirect effect in the extrapolation is not 
quite so negligible. In making evaluations based on 
the empirical constants of Benedict, it was found 
that gross differences in behavior between v 1 and 
v3 dependences were produced with vl or v3 
large, particularly as v2 was increased so as to be 
more than a small integer. This characteristic is at- 
tributed to the effect of the long range of the ex- 
trapolation with equations fitted to data at low 
quantum numbers only. A more uniform behavior 
has been obtained by a revised procedure for treat- 
ing the empirical vibrational energy. 



Energies for High Vibrational Levels 

In the last several years the method for estima- 
tion of high vibrational levels appears to be chang- 
ing, involving such new developments as are 
referred to as localized bond excitation and local 
mode description [10]. Special potential forms can 
be used for such calculations with constants con- 
verted [11] from empirical values found with a con- 
ventional valence bond system formulation and 
normal coordinate analysis [12]. Potential improve- 
ment based on direct comparison between com- 
puted and "observed" levels could be an ultimate 
objective. A hazard at the outset in this approach 
may be a sensitivity to the correctness of identifica- 
tion or assignment of spectroscopic data on which 
at least the original numerical constants are based. 
As to direct a priori quantum mechanical calcula- 
tion of levels for the molecule as a collection of 
nuclei and electrons, based on general physical 
constants, it appears that significant advances have 
been made on this intrinsically difficult endeavor. 



Somewhat approximate agreement with vibrational 
fundamentals has been obtained [13,14] but 
whether a similar quality of prediction could be 
achieved for higher vibrational levels may be in an 
area of pure speculation. A priori calculation ap- 
pears informative in regard to excited electronic 
states [15], in an energy domain beyond the range 
of the present treatment. 

Even if a local mode description will prove ulti- 
mately more reliable than the conventional ap- 
proach, it has appeared expedient to continue for 
the time being with the older formulation, for 
which the necessary parameters are at hand. It ap- 
pears plausible that newly and correctly calculated 
levels should on the average agree tolerably well 
with the old values of corresponding description. 
This is thought to be the usual situation for a group 
of "interacting levels" in a so-called resonance situ- 
ation. 

The vibrational constants used here are based on 
a formulation by the late Prof. W. S. Benedict [8], 
described by him as preliminary. His result can be 
shown as 

G(vl,v2,v3) = 3692.5965 vl + 1609.1113 v2 

+ 3803.6304v 3 -41.5442 V\ VI 

-13.4642 V2 V2- 48.0343 V3 V3 

-28.6309 V\ V2 

-164.2450 V\ V3 

-19.2960 V2 V3 

+0.0927 V\V\V\ 

-0.9003 V2 V2 VI 

+0.2690 V3 V3 V3 

-0.7760 Fl F2 F3 

+ 1.9316 V\ V\ V2 

+0.2325 VI V\ V3 

+ 1.0522 V\ ^2K2 

+ 1.1192 VI V1V3 

+ 1.5269 VI V3 V3 

-0.9318 V2 V3 V3, (48) 
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for levels with resonance shifts removed as indi- 
cated earlier. 

The present proposed innovation in regard to vi- 
brational energy is to suppress the long range ef- 
fects of Benedict's fitting on the basis that the fine 
details of fit while relevant in the region of fit in 
the low quantum number range (vi < 5) may still 
not be numerically reliable when extrapolated to 
large vi. For the various "small" quadratic and cu- 
bic terms, involving "v products" (=p), an extra- 
polation by replacement of "p" by p/[\+(p/a) k ] 
has been used with k=6 and with the parameter 
"a " chosen differently for terms quadratic and cu- 
bic in the v's. (65 versus 460) This causes these 
terms to become small in the approach to the disso- 
ciation region. 

For the estimation of vibrational levels in the re- 
gion of large v 1 and v3, the procedure adopted was 
to take the energy as given primarily by a quadratic 
jointly in vl and v3, much as in the case with a 
Morse potential in a diatomic molecule. Thus, in 
the case with v2=0, the form for this main part of 
the vibrational energy becomes 

GI(vl,0,v3)=i>l vl+v>3v3 -xllvl(vl-l) 

-x33v3(v3-l)-xl3vl v3. (49) 

The anharmonicities for this were chosen so as to 
agree with energies of dissociative reactions based 
on thermochemical data. For vl or v3 increasing 
singly with the other at zero, there is dissociation 
according to H 2 0=0+2H at about 76721 
cm~'=£>. For vl and v3 equal and advancing to- 
gether, dissociation is taken to be according to 
H 2 0=H + OH at about 41280 cm-'=Dm. 

In the cases of vl and v3 advancing singly, the 
energy is given as in 



G = i)v — x v (v — 1). 



(50) 



With i>1 = 3651.145 cm-' from G ( 1, 0,0) -G (0,0,0) 
and also with 03 = 3755.8651 cm"-' from G(0,0,1) 
— G(0,0,0), the corresponding "anharmonicity" 
constants follow from Birge-Sponer type relations 
as 



x = 2D - 1> - [(2D - 1>) 2 - i> 2 ] l/2 . 



(51) 



For D — 76721 cm' 1 , the long range estimates are 
^11=44.5045 and ^33 = 47.1274 cm"'. These are 
raised by about 0.0004 cm' 1 in covering small 
residual effects from the suppressed higher order 
constants at dissociation, which appears to be 
where vl or v3 singly reach a value of about 41. 



A somewhat similar procedure using Dm 
= 41280 cm' 1 to estimate xl3 for vl = v3 = v im- 
plies the relation 



G=i)v--xv ! 



(52) 



where v»=i>l +i>3+x ll+;t33 and x—x 1 1 -+-jc 33 
+xl3. The usual Birge-Sponer relations by 

x = i> 2 /4Dm or x 13 = i> 2 /4Dm -x 1 1 -x 33 (53) 

give x 13 =248.95779 kaysers but the partially sup- 
pressed residual contributions of other constants at 
this dissociation energy (near vl=v3 — 11) raise 
;tl3 to 251.3489 cm- 1 . 

A multiplier factor [1-0.0028 v2- 0.00013 v2 
(v2 — 1)] has been introduced for i>, vl+£i v3 to 
allow for a diminishing energy increment to disso- 
ciation as v2 advances upward above v2=0. All 
such adjustments are compensated for in the low 
quantum number range so as to preserve the behav- 
ior there according to the empirical data fit of 
Benedict [8], 



Thermodynamic Tables for H2O 

Two sets of tabular values have been included as 
prospective thermodynamic quantities for the ideal 
gas state of the light isotopic water molecule. 
These are here designated by their dates of compu- 
tation, which were 1982 and 1984. 

For the 1982 table, appearing here as table 2, 
parameter values used were influenced by results 
of computations for a rigid bender model, adjusted 
further for bond length increase by centrifugal 
stretching due to a rotational character of motion 
in the bending vibration. These included indica- 
tions as to the v2 dependence of the principal rota- 
tional constants, the extrapolation of vibrational 
energy to high v2 values, and the course of the 
ordinary centrifugal effects to high v2 and elevated 
temperatures, using a five parameter Pade formula- 
tion. 

The 1984 table, shown here as table 3, includes 
the innovations of the 1982 table, and a few others, 
also. In the approach to dissociation at high v 1 and 
v3, the behavior of Gv was taken as essentially 
quadratic in vl and v3, in resemblance to the 
known diatomic behavior with a Morse potential. 
Special functions were used to fade out the detailed 
higher order terms arising out of Benedict's Gv fit 
at low vibrational quantum numbers. For eight 
low-lying vibrational states, numerical values were 
inserted via the computer program for observed 
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Table 2, Thermodynamic quantities for light isotopic water (1982 version). 



T/K 

200 
300 
400 

500 
600 

700 

800 
900 
1000 

2000 
3000 
4000 

5000 
6000 
7000 

8000 
9000 
10000 



4.01111 
4.04064 
4. 12079 

4.23672 
4.36880 
4.50951 

4.65679 
4.80868 
4.96224 

6.17104 
6.78260 
7.13354 

7.38773 
7.60934 
7.79162 

7,90299 
7.92822 
7.87725 



5° 


G"-HS 


if 


RT 


21.09225 


17.11003 


22.72269 


18.72717 


23.89493 


19.87912 


24.82647 


20.77846 


25.61043 


21.52011 


26.29436 


22.15429 


26.90607 


22.71068 


27.46330 


23.20823 


27.97789 


23.65978 


31.83102 


26.85729 


34.46206 


28.97640 


36.46511 


30.60809 


38.08513 


31.94658 


39.45206 


33.08648 


40.63945 


34.08250 


41.68807 


34.96896 


42.62120 


35.76836 


43.45446 


36.49608 



h'-hS 

RT 



ficut 

Q 



3.98222 
3.99552 
4.01581 

4.04801 
4.09032 
4. 14007 

4.19539 
4.25507 
4.31810 

4.97373 
5.48566 
5.85702 

6.13856 
6.36559 
6.55695 

6.71911 
6.85284 
6.95838 



0. 
0. 
0. 

0. 
0. 
0. 

0. 

0. 
0. 

0. 
0. 
0. 

0. 
0. 
0. 

0. 
0. 
0. 



3 

6 

11 

17 
26 

37 

53 
72 

95 

635 
2366 
4519 

6524 
8080 
9487 

10818 
12075 
13255 



The final column gives the number of vibrational levels involved in the state sum. 



Table 3. Thermodynamic quantities for light isotopic water (1984 version). 



T/K 


/-i0 

R 


R 


C-nS 

RT 


H'-HS 
RT 


Gcut 
Q 


# 


200 


4.01111 


21.09218 


17.10996 


3.98222 


0. 


3 


300 


4.04065 


22.72262 


18.72710 


3.99552 


0. 


6 


400 


4.12080 


23.89486 


19.87905 


4.01582 


0. 


11 


500 


4.23676 


24.82641 


20.77839 


4.04802 


0. 


17 


600 


4.36895 


25.61038 


21.52004 


4.09034 


0. 


26 


700 


4.50993 


26.29435 


22.15423 


4.14012 


0. 


37 


800 


4.65779 


26.90615 


22.71063 


4.19552 


0. 


53 


900 


4.81075 


27.46355 


23.20820 


4.25535 


0. 


72 


1000 


4.96610 


27.97845 


23.65980 


4.31865 


0. 


95 


2000 


6.22473 


31.84805 


26.86066 


4.98740 


0. 


651 


3000 


6.83435 


34.50332 


28.98865 


5.51467 


3.290E-17 


2369 


4000 


7.13573 


36.51478 


30.62900 


5.88579 


3.328E-11 


4774 


5000 


7.34143 


38.12993 


31.97291 


6.15702 


2.054E-08 


7063 


6000 


7.52368 


39.48479 


33.11495 


6.36984 


4.151E-07 


8973 


7000 


7.67781 


40.65664 


34.11050 


6.54614 


9.326E-06 


10678 


8000 


7.77806 


41.68916 


34.99460 


6.69456 


3.320E-05 


11899 


9000 


7.80016 


42.60601 


35.79032 


6.81569 


2.384E-04 


12183 


10000 


7.76665 


43.42708 


36.51385 


6.91322 


4.060E-04 


12183 



The final column gives the number of vibrational levels involved in the final state sum. The next to the last column gives the 
fractional reduction in the state sum due to rotational cut-off, according to the "locus-asymptote" estimate used. 
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values for vibrational energy, principal rotational 
constants, and their five member Pade centrifugal 
parameters, based on reported spectroscopic data 
analyses using the Watson Hamiltonian formula- 
tion. A rotational cut-off approximation of a "lo- 
cus-asymptote" type was also introduced, but with 
little apparent effect up to 10000 K. 

It is natural to see the difference in values be- 
tween the two tables as relevant to their uncer- 
tainty. It is presumed that the disagreement in 
values should be attributed to effects in changes in 
level distribution, which may reflect the ad hoc 
modification of level description for the later table. 

It had appeared reasonable to maintain a favor- 
able view of progress in raising the number of con- 
stants based on the Watson Hamiltonian data from 
three to five. However, it is now recognized in ret- 
rospect that some basis for reserve exists. As used, 
the program for finding Pade constants was able to 
produce the five parameters as desired even when 
the Hamiltonian parameters were not complete to a 
corresponding extent. This might be termed a 
"spill-over" effect akin to the forming of product 
terms in a series development. The highest power 
of T directly included as a contribution to "T times 
the Wilson constant" may be obtained by taking 
the highest net power of J's in the Hamiltonian, 
dividing by 2, and subtracting 1. On this basis, the 
ground state and first excited vibrational state, 
(000) and (010), may be "complete" through the 
5th power. The states (020) and (030) show fitting 
in the 4th power, and the states (100), (001), (110) 
and (011) include only into the 3rd power. One 
may hope that a moving of the Pade process into 
the Hamiltonian will lead to a more uniform 
treatment [16]. 

As comment on our present use of a "pre- 
liminary" 1972 data formulation received from 
Professor Benedict [8], we accepted his view that 
his was better than that of Khachkuruzov [17], of 
1959. We note that a more recent vibrational 
energy formulation presumably of comparable 
quality was published in 1983 by Bykov, 
Makushkin and Ulenikov [18], and could in all 
probability provide a similar basis for a table of 
thermodynamic quantities. 

It appears that greater consideration should be 
given to recent work such as that by Child and 
Lawton [19] on local mode representations of 
vibrational states. However, at this time it is not 
clear how energies for the entire manifold of vi- 
brational states would be reliably and conveniently 
given for the calculation of thermodynamic func- 
tions on such a basis. 



Conclusion 

The objective in this study has been to obtain an 
improved extrapolation of the ideal gas table to 
higher temperatures. The procedure has made use 
of direct data, augmented with numerical estimates 
based on simple physical models. It is hoped that 
this might provide a realistic approach to better 
sum of state estimation. 

Although the models have involved some nu- 
merical choices that were not at all rigorous, the 
results may allow such comparisons as may lead to 
an informed appreciation of the problems remain- 
ing for the reduction of uncertainties. 
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APPENDIX 



A listing of the interactive computer program, in BASIC, for carry- 
ing out the preparation of the correction factor as a series in powers of 
temperature using the empirical constants of the Watson-type Hamilto- 
nian: 



10 PR8 



. : PRINT 



NOVEMBER 17, l-^&S 
'CENTRI4 1 OF H„ W. WOQ.LEV , MuV. 1^ 



-!IN7 ; PRIM 



£0 PRINT TSB< 8); "THE CENTRIFUGAL EFFECT OF MO._ECULP!R RDTH7ION ON IDS 

RL GOG PROPERTIES": PRtt 
30 PRINT "ENTER COMMENTS': INPUT C* : PRINT "ENTER NOME' ; INPUT 0»; iRINT 

"ENTflR DHTE" : INPU^ Df : PRINT LNrtR DDTP SUURCt" : INPUT H'F> 
40 PR# 1: PRINT TfiB< 2S > ; fl* : PRINT TOE ( 3a>;D*: PRINT TAB! 3>;L."1n PRINT 

Tft& < 5) ;S*: PRtt 
50 PRINT "ENTER RDTOTIONfiL CONSTANTS ft, B ,« C" : INPUT P!Z,F.X,CV 
60 DIM W* (30) , NT (35) , KF < 130) , NX 1 1 30) , NY (13*3! , NZ < 130) , H (6, 10) , CO ' 130 ), CF 

( 35 > , HH < 10 > , HL !1 i ', - < i 5 > 
7i2i DIM CNC40> , XX 140) , YY<40) , ZZ 140) 

33 FF i I ) = 1:FF<2> - 2:F-<3).^ £ : FF S 4 ) --^ 24:FF<5> - i*GUFFi&> -- 720 
90 PRINT "ENTER POWER UF T TO BE COVERED <>7>, OS l?" £3 3"": INPUT KN 
100 Fill = 1; FOB I = 2 "Q 15:F(I) = F(I 1 > « ( £ -» 1 3> / 2; >.t!T I 

110 SOTO 1 000 

200 REM iPROLIFFRRTICN PT BRANCH IMC 
ilB IK - 1:IL. -" lUCHi ■= 1:Il(KLi = II 
220 IK •' IK I 1: IF 1< > KL THEN 260 

230 IF ILriK) - ILiIK - i: THEN NC!IL? ■■ NC < J 'J > i~ i: CD TO ,22C? 
240 IU - IU + 1:NC<IU) - 1 
258 IF IK < KL THEN 22* 
360 IT = IU 
£712 FC = 1: FOR IU - 1 TO I T 1 F C - FC * FF INC < I U ; > : NEXT lUiCT -- FF'KL) / 

FC 
280 MX - MX + XXIII>:MY = MY + rY(Il):MZ = r.Z * 2 Z I I I i : MT ^ MX • «Y I- HZ 

- 3 - KL 
330 TR -^ - TR * CNCII) / KL : IF MT > KH THEN 310 
300 H<KL,MT) - H(KL, MT) *- TR * F!MX) » F!MY> » F<M7) / P1Z !MZ - .) • 3 

X ■ <M> -1) / CY ' (MY - 1) * CT 
310 RETURN 

400 REM ;RL ! URNING IQ fi PlSC£ UF BRANCHING 
410 MX - MX - XXSII):MY - MY - YY(II):MZ - MZ - ZZ(III:TS = - TR * KL / 

CN I I I ) : RETURN 
1000 REM 
1010 W*(l> - "DELJ" :W<M2> - "DELJR" :W*<3> "• "DELK" : W* (4 > - 'SDLJ";W*<3> - 

"SDUi" :U»iE) = "HJ"!U*(7! - "HJK" : W* CS> - "HKJ" 
1020 W*(9) - »HK":W*(ie> = "SHJ" :W*(11) ~ '■ SHJK ■' : W* { IE) - "SHU" : W* 1 13) - 

"CL-" ;W*l 14) = "CLK":H*(15I -- "CLKK J " : W* < 16) = "CLJ"." 
S030 H*<17; -= "CLJJK" :W*<18> = "SL J " : W* ( I 9) ~ "SLK" :Ui«en - "SiHj":U»« 

1) = "SLJK" :W*<2S) "- "CPK" 
1040 W*<a3) = "CPKKJ" : W*(P4> = "OPKJ " : iJ$lK ) ■* "SPK " ; US* <2E) - "Z.K~":W*1F 

7> = "Z10P2":k«!28) = " Z 1 4 " : W* 1 ETJ) ■■- " Z 1 1 " ; W* ( 36" ) ■■- "ZlB" 
1050 DfiTfi 1, S, - 1 , E 1 , 0, iZt, -1 , 0, E, 0, -1, 0, 0, C, t, 1, 1,*, -^, 1,?, 1, £,&,!,! 
10E0 DfiTfi E, 3, -1, 1, 0, 1 , -1 , 0, 1, 1, - l, 0, 0, E 
1070 DOTO 3,1,-1,0,0,2 

1060 DRTFR 4, 4, -£, 3, 0, 0, 3, 0, .^, 0, -F, 1 , 0, 1, E, 0, 1 , 1 
1090 DftTfl S, 3, -3, 1, 0, 1, a, 0, 1, 1 
1 I 00 DATft 6, 10, 1,3, 0,0, 1,0, J, 0, 1,0, 0, 3, 3, 3,1,0, 3, 3, 0, 1,3, 1,3,0, 3, 0,3:, 1 

, 3, 1, 0,3, 3,0, 1, a, 6, 1, 1, 1 
1110 DfiTP, 7,b, 1,^,0, 1, 1,0, E J , 1, 1,0, 0, 3, 3, 1, 1, 1,3, 1,0, 3, 3, 0, 1,3 
1130 DP1TR 8,3,1,1,0,3,1,0,1,3,1,0,0,3 



1130 DfiTP, 9,1,1,0,0,3 

1140 DRTO 10,8, 3, 3, 0, 0, 3% 1, 3, 0, 3, 1, 0,3 , 3, 3, 1 , 0, 4, p, 0, I , - .'■■., 0, 3, 

, -3, 0, 1,3, -4, 0, 3, 1 
11LJ0 DATA 11,4,3,2,0,1,-3,0,3,1,3,1,0,3, 3,0,1,3 
1160 DOTO 12, 3,3, 1, 0. 3, -S, 0, 1, 2 
1 170 DFiTA 13, 15, 1, 4, 0, 0, 1, 0, 4, 0, 1, 0, 0, 4, 4, 3, .1, 0, 4, 1, 3, 0, 4, 3, 0, 1 , 4, 0, 3, 

1, 4, 1, 0, 3, 4,0, 1, 3, 6, 2, 2, 0, F, 3, 0, 2, E, 0, 3, 2, 13, 2, 1, 1, 12, 1,3, 1, 12, 1, 1 , 



1,0,0,4 
, 3, 1 , 1 , 0, 3, 1 , 0, 1 , 



1, 1, 3, 2, 1,0, . 
: , 2 , 1 , 1 , 3 , 2 , 



, 0, 1, 3, 
- fa, 0, 3, 



-2, 0, [ , ' 
2, 1, 2, 1, 



3, 1, 2, I , 3,0, 2 
I , 0, H, 3 , 3, 0, C, 



1, 0, : 



-4, 0, 



, 0, 1,4 



-2, 



TERMS: 
11);") 



1180 DOTO 

1190 DP)TH 

1300 DOTO it,S, 1, 2, 0, 2, 1, 0, 2, 2, 1, 0, 0, 

1210 Df-lTA 17,10,1,3,0,1,1,0,3,1,1,0,0,4 

, 2, 3, 1, 0, 3, 3, 0, I , 3, E, 1, 1, 3 
1220 DOTA IS, 13, 3, 4, 0, 0^ -2, 0, 4, 0, Z 1 , 1,2,3 
3, 0, 1, -E,, 0, 3, 1, &, 2, 1, 1, -€,, 1, 3, 1, 6, 2, 
1230 DOTO IS, 2, 2, 1,0, 3, -3,0, 1,3 

1240 DP1TO 20, 4, £, 2, 0, 2, -2, 0, 2, P, 3, i , a, I 

1350 DOTP) 21, S, 2, 3, 0, 1, 2, 2, 1, 1, ■», 2, 

, 1, 3, -3,0, 3, 1 
IE'20 DftTft 22,1,1,0,0,5 

1270 DfiTA 23, 3, 1, 1, 0, 4, 1, 0, 1, 4, 1, 0, 0, S 

12S0 DOTO 24, £, 1, 2, 0, 3, 1, 0, 3, 3, 1, 0. 0, 5, P, 1, 1, 3, 1-, 1, 0, < 

12-30 DPTP 3^,2,2,1,0,4,-2,0,1,4 

1300 DOTO 26,1,1,0,0,6 

1310 DP)Tfi 27,3,1,1,0, S, 1, 0, [, 5, 1, 0, 0, F, 

1320 DRTft 28,1,1,0,0,7 

1330 DOTO 2S, 1,1,0,0,S 

1340 DOTft 30,1,1,0,0,9 

1350 DRTft 0,0,0,0,0,0,0,0,0,0, 
13&0 1=0: SPtbU= 200 

1370 REfiD NN, N9: IF NN = T-HEN I4E0 
1380 PRINT : PRINT 

1390 FOR Kl = 1 TO N3:l = I J- 1 

1400 REOD KFtI ) , NX< I) , NY < I! , NZ II) : PRINT ■'»" 
;KF<!);", (NX,NY,NZ): ( ' ;NX I I > i " , " ; NY < 1 ) 
NN) - N9:LX = NN 
1410 60TU 1370 
1420 KN = I : SPFFD= PS5 
1430 FOR N = I TO LX: PRINT "ENTER COEFFICIENT ' C= ( " 

INPUT CF IN) : NEXT N 
1440 PR# 1: PRINT " CONSTPNTS IN THE RQTfll IONfiL HOM1LTONIOM" : FOR U *• I 1 
LX: PRINT CF(N),W*(N): NEXT N: PRINT " R* " ; fiZ ; " . B=";F<X;" r=-"-cv- 
PR# ' 

1450 1=0: FOR N = 1 TO LX:NB - NT(N) 

1460 FOR K = 1 TO N8.-I =. I + l:CO<t) « KF ! I > » CF<N); NEXT K: NEXT N 
1470 FOR J = 1 TO 40:CN<j; • O: NEXT J 

1480 JM = 1:CN(1) = C0(1):XX<1> = NX(H:YY!1) = NYC):ZZ(1) - :1Z '. 1 ) 
1490 FOR I = 2 TO KN:KK = 0: IF TDdi - THEN 1 5E0 
1500 FOR J = 1 TO JM 
1510 IF <XX(J) - NX(I>) ■ 2 + <VY(J) - NY ( I ) ) ■■ 2 * <ZZ<J> MZ ! I ) ) 

- THEN KK - 1; GOTO 1530 
15P0 GOTO 1540 

1530 CN<JI - CN<J) t- CO(II:J = JM 
1540 NEXT J 
1550 IF KK * THEN JM = JM + 1:CN(JM! = CO C I ) : X X ( JM ) - NX<I):YY(JM> = 

NY! 1 ) : ZZ (JM) = NZ ( I) 
1550 REM 
1570 NEXT I 
1580 FOR J = 1 TD JM 
15^0 PRINT "*";J;", <XX,YY,ZZ): ( " ; XX ( J) ; ■■ , » ; YY ( J) : ■■ , » ; ZZ < J > : " ) " ; " , COF 

F: ";CM<J> 
1600 NEXT J:KN =- JM 
1610 FOR MT = 1 TO 10: FOR KL = 1 TO 6:H(KL,MT) = 0: NEXT KLsHH(MI) - 

: NEXT MT 
1620 MX = 1 :MY = 1 : MZ * 1 : TR = 1 : MT = 

1630 FOR 10 = 1 TO KNiKL = 1:11 = ifl. IF ctM ( 1 1 ) = THEN l-#3B 
1640 PRINT "I1=";I1;", TR=";TR 



:f«;", COEF: 
: NLK1 K I :,\>T 



' ;Wt 



16SB 
I860 
16721 
1600 

i&ga 

172113 
17121 
1720 
17351 
174B 
1750 
17613 
1770 
1780 
1718 
1800 
1810 

iosa 

1830 
1848 
1850 
18t,0 
1870 

ieaa 

1890 

1900 

1910 
1920 
1930 
1940 
1950 
I960 
19 70 



1990 
£0013 
£0 1 



2030 
=2040 



IF MT ! KH THEN Ifl = KN : GOTO 1<330 
11 = II: BOSUB £00: IF KH < KL + 1 THEN 1320 

FOR IB = IS TO KN:KL = £:IE «= IE: I- CN<I£) =-- a THEN 1910 

IF MT > KH THEN IB = KN: 6DTD 1910 
II = I£: GOSUB £00: IF KH < KL - 1 THEN 190a 

FOR IC = IB TO KN:KL = 3:13 = IC: IF CNCI3) = B THEN 1B90 

IF MT > KH THEN IC = KN: EOTO 1S90 
IT = 115: GQRUB £00: IF KH { KL + 1 THEN 1880 

FOR ID = IC TQ KN:KL > 4:14 = ID: IF CN(I4) = THEN 1 07lZl 

IF MT > KH THEN ID = KN : GOTO 1870 
IT = I4i GOSUB £00! IF KH < KL » 1 THEN 1S60 

FDR IE = ID TD KM : KL - 5: IS = IE: IF CNII5) = THEN 1 B50 

IF MT ) KH THEN IE -- KN : GOTO 1850 
II - 15: GOSUB £00: IF KH < KL t 1 THEN 1S40 

FOR IG = IE TU KN:KL = 6:16 = IG: IF CM (16) '= 8 THEN 1B30 

IF MT > KH THEN IG J- KN: GOTO iS3« 

I 1 = 16: GOSUB £00 

II = IS:KL = 6: 60SUB 400 

NEXT IG: Ml • MX + MY + Bi - e - KL 
II ■= I3:KL = 5: GOSUB 4«iei 

NEXT Ib:«l -' BX t <IT i-W - e - KL 
[I = 14:KL = 4: GOSUB 400 

NEXT IC;!tI - MX + MV + MZ - E - KL 
II ■= I3:KL - 3: GOSUB 400 

NEXT IC-.MT - MX + MY + MZ - £ - KL 
II = l£:KL = £: GOSUB 400 

NEXT IB:MT = *IX + MY + MZ - £ - KL 
II = I1:KL = 1: GOSUB 400 

NEXT Ifl 

FOR MT = 1 TO KH: FOR KL - 1 TO KH 
HH(MT) = HH1MT! + H(KL.MT) 

NEXT KL: NEXT MT 

PR* 1: PRINT : PRINT "COEFFICIENTS OF PROPORTIONAL ADDITIONS TO O, 

BY POWERS OF T/C2" 

PRINT "HH: 1 : " jHH(l) ; ";£: " ;HH<£> ; " ;3: ";HH!3> ; " ; 4 : " ; HH (4) ; " j 5: " ; HH ( 5 
) ; " : 6 : " HH ( 6 > 

GOSUB 5000 

PRINT "ADDITIONS TO LOG 0, BY POWFRS OF T/C£" 

PRINT "COEFFICIENTS:": FOR N = 1 TO KH : PRINT "POUER=" ; N; ■' , COEF, 
=" ;HL (N) : NEXT N 

C£ = 1.43S7B6: FOR N =» 1 TO KH;fi<M> ■= HL<N) / C£ ' N: PRINT "COEF. 
OF T";N;" a";n(Nla NEXT M 

IF KH < S THEN £090 
Bl = Pl<£> / 0<1> ;B£ - Ft (3) / fi < 1 > : B3 = 0<4> / O < 1 ) : B4 = 0<») / tt < 1 > 



v S + £ * HHC3) * HH<5) * HH<1) 
HHI£) - £ ♦ HH<1) - £ * 3 / £ i 



4) * HH(1) 

HH ( 1 ) " 3 

1 ) ' 6 / 6 
51 30 RETURN 

6000 PR* 1: PRINT TSB ( 039) ;D»: PRINT 
6010 PRINT CHRt <9);"60P"; 
6020 LIST : PRINT TAB ( 039) ;D*: PR» 0: END 

NOVEMBER 17, 



HH (£> 
HH ( £ > * 



: 3 / 3 - HHC3) * 
HH ( 1 ) ' 4 - HH < 



CHR* <3> ; "80N" 



3:D1 = (Bl 



B£ 



*i350 DN =-- B£ ■ £ Bl * 

D£ » B4) / DN 
£0£0 Gl - Dl -i- B1:G£ - D£ + Dl * Bl + B£ 
£070 61 = Gt * 0(1):6£ - G£ * 0(1) 
£080 PRINT : PRINT "D LN Q = < " ; ft < 1 > ; " * T + (' 

) * T ' 3) / (i + (";D1;") * T + C *' ; D£ ; " ) ' 
£030 PR# 0: END 
5000 REM LOO 
5010 HLC1 ) = HH(l) 
50£0 IF KH < £ THEN 5150 
5030 HLCS) = HH(L') - HH ( 1 ) "- 2 / 2 
5040 IF KH < 3 THEN 5150 
5050 HI. (3) = HH<3) - HH < 1 ) 
5060 IF KH (4 THEN 5150 
5070 HL(4! = HH(4) - HH(E) 

HH < 1 ) '• 4 / 4 
50S0 IF KH < 5 THEN 5150 
5090 HL<5) = HH(5) - HH(2> 

HH(1) " £ * HH<3) - HH<£> 
5100 IF KH < 6 THEN 5150 
5110 HL<6) = HH<6) - HH(3! -■ 2 / 



BE) / DIM: DC - <B3 



31 ; " > 



HH(£ 


) + 


HH(1 ) 




3/3 




2 - 


HH < 1 ) 


* 


HH(3> 


HH<3 


;> - 


HH(l) 


* 


HHI4) 


:> * 


HHS1) ■' 3 


+ 


HHU) 



HH<5) <- HH( 



